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1 Introduction 



The present experimental data on lepton mixing |TJ |2] are compatible with maximal 
atmospheric neutrino mixing. One possible theoretical explanation (for recent reviews 
see jH]) for that maximal mixing, if it indeed occurs, is /x-r interchange symmetry [HE]. 
In the basis where the charged-lepton mass matrix is diagonal, we denote the effective 
light-neutrino Majorana mass matrix by Ai u . The diagonalization of M. v is performed as 



U T M v U = diag (mi, m 2 , m 3 ) 



where mi, 2 ,3 are the neutrino masses and U is the (unitary) lepton mixing matrix. If A4 t 
is fjL-r symmetric, i.e. if M. v = with 



M< s > 




(2) 



then atmospheric neutrino mixing is maximal, i.e. \U^\ = \U T ^\. Indeed, the matrix (0) 
has six physical parameters, since only the phases of zw* and of y 2 x*z* are physical. 
The existence of only six physical parameters in a theory for nine observables — the three 
neutrino masses m 12i 3, the three lepton mixing angles 6*23,13,12, one Dirac phase 5, and 
two Majorana phases and Q — leads to three predictions. 1 In this case the predictions 
are the following: {U^l = \U T z\ (#23 = 7r/4), U e % = (#13 = 0), and the Dirac phase 6 is 
meaningless since U has one vanishing matrix element. 

The fi-r interchange symmetry, represented by the matrix 



T 




(3) 



is a Z2 symmetry, and as such it has two possible eigenvalues: +1 and —1. The matrix 
j\/f(S) corresponds to the eigenvalue +1: TM^T = +M^ S \ But there is also the possibility 
that M u is /i-r antisymmetric 0, i.e. that M v = M^ AS \ where TM^T = -M^ or 




M (AS) 

This possibility seems attractive since M^ AS ) may be diagonalized as 

yT M (AS) y = diag ^ _ ik ^ ^ 

where k - 

( 2y* 2y* -2z 



(4) 



(5) 



2 I y 1 2 + I z 1 2 and the unitary matrix V is 



V 2k 



z* +ik z* — ik 2y 
\ z* - ik z* + ik 2y 



(6) 



1 We defer the precise definition of the mixing angles, Dirac phase, and Majorana phases to section 3. 
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Thus, M. v = M( AS ) leads to two degenerate massive neutrinos, one massless neutrino, 
maximal atmospheric neutrino mixing, and maximal solar neutrino mixing (\U e i\ = \U e 2\ 
or 612 = 7r/4). These predictions are not so far from reality, if the neutrino mass spectrum 
is inverted. However, the exact equality Ai v = M^ AS ^ can be excluded, since we know 
that all three neutrino masses are non-degenerate and that solar neutrino mixing is large 
but not maximal. 

In this paper we consider various possible perturbations to /i-r antisymmetry. More 
precisely, our starting point is 

M -l = M (AS) + M (P) j (7) 

where is a /i-r symmetric matrix, a "perturbation" to M^ AS \ Having A4.~ l instead 
of M. v in equation (J2J) is quite advantageous for model building; this will become evident 
in section 2, where we construct a class of models based on equation (J2J). Physically, 
exchanging M.~ x with M. u simply corresponds to complex-conjugating U and trading 
the neutrino masses by their inverses. In our class of models, the perturbation is 
mostly free, since it depends on the assumed scalar content of each model. We proceed 
in the following sections to analyze in detail the consequences of some perturbations. We 
show analytically in sections 3-5 that the simplest possible perturbations lead to very 
constrained neutrino mass matrices, which contradict the experimental data. We proceed 
in section 6 to analyze numerically some more complicated perturbations wherein A4 U 
has six physical parameters. We summarize our findings in section 7. 



2 A class of models 

Consider a SU(2) x £7(1) gauge theory with three scalar SU{2) doublets 4>j {j — 1,2,3). 
Besides the left-handed lepton SU(2) doublets D a and the right-handed charged-lepton 
577(2) singlets a R (a = e, /i, r), we introduce right-handed neutrinos v aR . We assume the 
validity of the lepton-number U(l) symmetries L a for all dimension-4 terms in the La- 
grangian; dimension-3 terms are allowed to softly break the L a . We assume the existence 
of a symmetry under which 0i, e R , and all the v aR change sign, while all other 

fields remain invariant. The Yukawa Lagrangian of the scalar doublets 2 is then 

A> = - (yiD e v eR + y 2 D fl u llR + y' 2 D T v T Rj iT 2 <t>X - y^D e e R (j)i 

- (vaD^r + y' 4 D T TR) 4> 2 - (ysDnfiR + y' 5 D T r R ) <f) 3 + H.c. (8) 

In the /i-r symmetric model built by two of us a few years ago [7j, a /i-r interchange 

(tr) 

symmetry 7L 2 was introduced, under which «-> D T1 fi R «-> t r , <-> v tR , and 

03 — > —03, such that after two applications of 7L 2 all fields transform into themselves. 
Instead, in the present model we employ a symmetry Z 4 under which 

D e -> iD e , Df, -> iD T , D T -> 

eR -> ie R , fi R -> ir R , t r -> ifi R , , g , 

VeR -> WeR, ~^ W T R, V tR -> Zf^, 

03 -» -03- 

2 We shall soon introduce some scalar SU(2) singlets into the theory, which will have further Yukawa 
interactions. 
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ftrl 

Notice that all lepton fields change sign after two applications of this Z4. Just as 7L 2 
in [7j, the present Z 4 enforces 2/2 = 2/2,2/4 = 2/4, and y' 5 = —2/5. The neutrino Dirac mass 
matrix Md and the charged-lepton mass matrix Mi are then 

M D = diag(a, b, b) , (10) 
M £ = diag (e^ e m e , e^m M , e^m r ) , (11) 

respectively, where a = y*Vi, b = JJ2V1, m e = \ y 3 vi |, = 12/4^2 + 2/5^31, an d J7i T = 
(2/4^2 — 2/5^3 1 , the Vj being the vacuum expectation values (VEVs) of the lower components 
of the 4>j. 

Now consider the Majorana mass terms of the right-handed neutrinos, 

1 I UeR \ 

^Majorana = £ { V eR, rfn> "tr) C' 1 M* r \ V^r + H.C. (12) 

V V rR J 

They have dimension 3 and are, therefore, allowed to break the lepton-number symmetries 
L a . Because of the symmetry Z 4 , which is not allowed to be broken softly, one has 

M R = M (AS) . (13) 

The seesaw mechanism [BJ El prescribes 

M v = -M T D M R l M D , (14) 

hence 

M~ l = -Mv X M R M T D \ (15) 

It follows from equations (fTU)) and (fTB^) that M.^ 1 is \i-t antisymmetric, just as Mr. 

The model as it now stands is not realistic for many reasons, in particular because Mr 
in equation (fT3*|) is not invertible and therefore the seesaw formula (fT4*j) cannot apply. One 
may correct for this by adding to the model one or more scalar SU(2) x U(l) invariants 
with non-vanishing lepton numbers. The simplest possibilities are the following: 

1. A complex scalar \ ee with L e = —2 and = L T = 0, which changes sign under Z4. 
Its Yukawa coupling Xee^eR^~ lu eR + H.c. generates a non-vanishing (M#) ee upon 
Xee acquiring a VEV. One then has 



case 1 : M. v 




(16) 



2. A complex scalar x^r with L e = and L M = L T = —1, which changes sign under Z4. 
Its Yukawa coupling X^t^rC^^tR + H.c. generates a non-vanishing (M R )^ T upon 
X^t acquiring a VEV. One then has 



case 2 : M^ 1 = \ y z w | . (17) 
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3. Two complex scalars Xe/i and Xer, Xe/i having L 
has L P = L T = — 1 and L 



= — 1 and L T = while Xer 
... - 0. Under Z4, Xe^ ^ — Xer- The Yukawa couplings 

V Ir C ~ X {Xe^nR + XerVrR) + H.C. lead tO (M R ) e/I ^ - (M R ) eT Upon Xe M and Xer 

acquiring VEVs. One thus obtains 



case 3 : M,, 1 = 



( 

y 



z 








(18) 



4. Two complex scalars Xnn and Xtt, having L e = L T = and = —2 while Xtt 
has L e = = and L T = —2. Under Z4, Xnn ^ ~Xtt- The Yukawa couplings 
Xw u Ir C ~ 1u hR + Xtt^ r C~ x v t r + H.c. lead to (M fl ) ^ - (M r ) tt . One then has 



case 4 : .M, 




(19) 



One may put together any two of the possibilities 1-4 and obtain a neutrino mass matrix 
with more degrees of freedom — and correspondingly less predictive power. One thus 
obtains six more possibilities: 



case 5 



case 6 



case 7 



case 



case 9 



case 10 



Ml 



M 



-1 



Mi 



Ml 



Ml 



Ml 




(20) 



(21) 



(22) 



(23) 



(24) 



(25) 



It is the purpose of the rest of this paper to study, both analytically and numerically, 
whether cases 1-10 are viable or not, and to find out their predictions. 

The neutrino mass matrix M u may be rephased: M u — > XM. V X, where X is a 
diagonal unitary matrix. Because of this possibility, in general only three phases in the 
symmetric M. v are physically meaningful. Using this rephasing freedom, one finds that 
cases 1-4 have four physical parameters each, case 10 has five physical parameters, and 
the remaining cases 5-9 have six physical parameters each. 
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3 Cases 3, 4, and 10 are not viable 

These three cases contradict experiment since in all three of them the conditions 

(M- 1 ) = 0, (26) 



(X 1 ) = (27) 

simultaneously hold [Ej. Using 

M' 1 = C/diag (mf 1 , m^ 1 , mj 1 ) U T ', (28) 
we find that the system of the two conditions and (|27jl leads to 

m 1 u el u r2 u* r3 - u e2 u^u; 3 



m 3 UaUrtUh-UaU&U*^ 
m 2 U e iU T2 U* 3 - U e2 U l j,iU£ 3 



(29) 
(30) 



m 3 u e3 u Tl u* 2 - u el u^u; 2 

We use the standard parametrization 

U = P F UP M , (31) 

where 

(Cl3 c 12 C13S12 Si3e 2 \ 

— C23S 12 — S 23 S V3 Ci 2 e lS C23C12 — S23 s 13 s 12e 1 ' 5 S23C13 , (32) 

—s 23 S\ 2 + c 23 S\ 3 C\ 2 e l S23C12 + C23S 13 si2e i — C23C13 / 

P F = diag(>\ e**", e^) , (33) 

P M = diag(e* e / 2 , l,e^ 2 ). (34) 

In equation (|32j) . s*. = sin^ and = cos^ for = 23,13,12. The phases d a in 
equation (J3~3~j) are unobservable. Using the parametrization above, equations and (JSUj) 
translate into 

m l gi(e-fl) ~ l C | 2 ) [ C 23-523 ("5 2 2 _ C 12) + C 12-Sl2 (-s| 3 - C23) e] 

m 3 c 23 S23S 2 2 (l-2 |e| 2 ) + c 2 3S23C 2 2 e* 2 + c 12 s 12 (c 2 23 - s 2 23 ) |e| 2 e* ' 

m 2 „_<n 



1 - |e| 2 ) [c 23 s 2 3 (s 2 2 - c 2 2 ) + C12S12 (sla - cf 3 ) e 



m 3 



e " 7 — " 2\ — ' ( 36 ) 

-C2 3 s 2 3C 2 2 (1 - 2 |e| J - c 2 3S23sf 2 e* 2 + c 12 s 12 (c| 3 - s 2 23 ) \e\ e* 



where we have defined e = Si3e* 5 . Experimentally [I], both |e| and \c 23 — s 23 \ are at most 
of order 0.2; neglecting terms of order e 2 , (c 23 — s 23 ) 2 , and e (c 2 3 — S23) in equations (J33j) 
and (|3l))). one obtains 

^ » 1-% (37) 
m 3 S12 



(38) 

m 3 cf 2 



y»2 _ 1 s 12 

r,2 
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Using the experimental result PQ tan 2 # 12 ~ 0.43, one obtains 



ml 



1.9. 



rrin — m 



(39) 



Although the number in the right-hand side of equation ()39|) is rather sensitive to the 
precise value of tan 2 # 12 , its order of magnitude does not change. It is in manifest contra- 
diction with experiment, which gives 



2 2 

ml — rn{ 



0.037. 



(40) 



Thus, cases 3, 4, and 10 are excluded. 



4 Case 1 is not viable 

Define the orthogonal matrix 

Ri 



I 1 

r r 
\ —r r 

Then, with the Ai' 1 of case 1, in equation (fTHjl. one has 



where r = — =. 

\/2 



(41) 



RjM^Ri 



(42) 



Z v^y 
v^y z 
V ^ 

This is a mass matrix of the Fritzsch type, which may be diagonalized in a standard 
fashion [TT]. The three phases in the matrix (|42jl may be discarded by using a diagonal 
unitary matrix Pi: 



P 1 R^M- 1 R 1 P 1 



I \x\ 

V o 



V2\y\ 
\z 
\z\ 



(43) 



This matrix has eigenvalues Ai, — A 2 , and A 3 , with Ai > A 2 > A 3 > 0. If we define the 
diagonal unitary matrix 

P 2 = diag(M, 1), (44) 
then there is an orthogonal matrix i? 2 such that 

P 2 RlP 1 R T 1 M- 1 R l P 1 R 2 P 2 = diag (A 3 , A 2 , Ai) . (45) 
The standard diagonalization of the Fritzsch mass matrix yields 



Mil 
(^2)12 
(^2)13 



A 3 (Ai + A 3 ) (A 2 - A 3 ) 
\ (A1-A3) (A 2 + A 3 ) (Ai-Aa + As)' 

A2 (Ai — A 2 ) (A 2 — A3) 

\ (Ai + A 2 )(A 2 + A 3 )(A 1 -A 2 + A 3 )' 



Ai (Ai - A 2 ) (Ai + A 3 ) 
\| (Ai + A 2 ) (Ai - A 3 ) (A!-A 2 + A 3 )' 



(46) 
(47) 
(48) 



Since 



in 



-i 



rn 



-i 



WM^U* = diag(m^, 

we see that there are two possibilities: 

1. If the neutrino mass spectrum is normal, nil < m 2 < m 3 , then Xj 
j = 1, 2, 3 and 



u e2 


2 


(-^2)12 






.(^2)13. 



m 9 



m 2 1 



m- 



m 



1 (mj 



+ m- 



m 2 1 + m 3 1 



(49) 



for 



(50) 



2. If the neutrino mass spectrum is inverted, m 3 < m% < m 2 , then Ai 



m 1 , and A3 



m 2 . One then has 

2 



u e2 


2 








.(•^2)12. 



m. 



1 fa 



- 1 ) K 1 



+ m 



-1 



m 1 1 (m 3 1 — mj x ) fm 3 



(51) 



m 



The inverted case is readily excluded. With m 2 — m\ = 8.1 x 10~ 5 eV 2 and \m\ 
2.2 x 10 _3 eV 2 , one finds that equation (|5Tj) yields |t/ e 2/t4i| 2 > 0.98 for all m 1; which 
is far above the best-fit value \U e2 /U e i\ 2 = 0.43. As for equation (|50|). with the above 
values for the mass-squared differences, it always yields \U e 2/U e i\ 2 < 0.28. This value 
translates into s\ 2 < 0.22. Reporting here for further use the 3a values derived in jT], 



7.2 x 10~ 5 eV^ <m\-m\< 9.1 x 10~° eV z 



1.4 x 10' 



; eV 2 < 



n? 



< 3.3 x lO^eV' 



0.23 < sin 2 # 12 < 0.38, 
0.34 < sin 2 # 23 < 0.66, 
sin 2 13 < 0.047, 



(52) 
(53) 
(54) 
(55) 
(56) 



we see that 0.22 is just outside the 3a range of sin 2 #12. However, the upper bound on 
\U e2 /U e i\ 2 of equation (JSTJj) is sensitive to the values of the mass-squared differences. 
Taking the lower limit m\ 



m 



7.2 x 10 5 eV from inequality (|52|) and the upper limit 
m\ = 3.3 x 10~ 3 eV 2 from inequality (J53J), the upper bound on \U e2 /C/ el | 2 is 0.35, or 



s 2 2 < 0.26, which is inside the 3a range of inequality (|54p. In that case, a simultaneous 



consideration of s? 



13 



[(i?2)n] 2 leads to s 2 3 



0.29, for the m\ where the maximum of 



equation ()50|) is reached. In view of the inequality ([56)1 . this value is clearly much too 
large, hence case 1 is excluded. 



5 Case 2 is not viable 

Case 2 is defined by 



K 1 ), 



0. 
0. 
0. 



(57) 
(58) 
(59) 
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In the following, we shall use the quantities 

Hi = mxe" 40 , /i 2 = m 2 , H3 = m 3 e~ tn . (60) 
The matrix elements {M-Z l ) a p are then given as 

{M-'U = E (61) 

3=1 y-j 

with U in equation We shall moreover use 

P = #r- (62) 

e = s^exp and v = cos 2^23- We also introduce the small parameter A = 0.22 which 
is useful since, experimentally, both |e| and \v\ have upper bounds of order A. Moreover, 
(m 2 — m\) j \m\ — m\\ ~ 0.037 is of order A 2 . 
Equation (joTj) yields 

(Uel) 2 (Ue2) 2 (Uesf , , 

^ '— + ^ '— + ^ '— = 0, (63) 

A*l A*2 ^3 

or 

^ = - r ^f^ + ^V 1 - (64) 
1 - | e | \fii Hi) 

Taking into account that m 2 > mi and ci 2 > si 2 , equation (JH%|) leads to the inequality 

m7 1 1^ r 2 - s 2 ' ( } 

mi i — |e| c 12 s 12 

hence the neutrino mass spectrum in case 2 must be inverted. The quantity c 2 2 — s 2 2 is 
smaller than 1: with the best-fit value s\ 2 = 0.30, it is 0.40; at the upper edge of the 3a 
range for s 2 2 , it is only 0.24, i.e. of order A. Thus, it is appropriate to admit m 3 /rrii to 
be of order A, instead of A 2 . As for rri3/m2, it is of the same order as 777.3/7711, since 

^ = , m3 (66) 
m 2 ym 2 + Am| 

and m\ is of the order of the atmospheric mass-squared difference when m 3 /mi is very 



small. We obtain a picture in which 777,1 ~ m 2 »s y Am 2 tm , while m 3 /mi < 0.2 is small. 
Equation (J55J) gives, after using equation 

s 23 - e i/3 c 23 + C (c 23 + e i/3 s 23 ) = 0, (67) 

where 

C = e -5 ; — c i2-si2- (68) 

Ci 2 /X 2 + s( 2 Hl 



Next we address equation (JH9*|) . Defining 

c 23 e *' 

cos /? + sin /3, (69) 



21 — c 23 e ^ ~l~ s 23 e 
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we find that equation may be written as 

az + bz* + c = 0, (70) 



where 



a = l-2|e|\ (71) 
b = -e* 2 ^ + ^1 , (72) 

Ci2/^2 + sf 2 /ii 



c = 2i|e|Vl-^ 2 Csin/?, (73) 

where we have once again used equation ()64j) and the definition (|68|). The solution of 
equation (J7U|) for z is 

be* — a*c ,_ 

* = | ,2 |,,2 - (74) 

|a| — |o| 

An upper limit on the absolute value of ( is 

IAI / II ^2+^1 

ICI < e-i 2 C 12 S 12 

Ci 2 m 2 - sf 2 mi 

2 (m 2 + mi) 

C12S12 



(c 2 2 - S12) ( m i + m 2) + m 2 - m 1 
< \e\ tan 20i 2 . (75) 

The 3cr limits |e| < 0.22 and tan2#i 2 < 4.04, cf. inequalities (|H4*|) and (JSHJ), lead to 
ICI < 0.89; at the 2a level — see |T] — this bound is already down to |£| < 0.50. 

Let us now use simple approximations to show that case 2 is not viable. Since a ~ 1 
and |6| < A 2 , we conclude from equation (f7i|) that 2 ~ — c: 



cos /? + sin j3 ~ — 2z \e\ 2 \/l — v 2 £ sin /3. (76) 
The right-hand side of equation (|76|) being at most of order A 2 , one must have 

e 4/3 ~r]i(r} = ±1), 1/ ~ 0, (77) 
up to corrections of order A 2 . Therefore, 

c 2 3 + e^23-^, s^-e^^ 1 —^. (78) 

Introducing the approximations (|78|) into equation (JB7J) . one finds 

C - = ( 79 ) 

1 + 27] 

i.e. C = ±z. But we know that at 3a level |£| is at most 0.89. Hence case 2 is ruled out. 
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6 Cases 5—9 are viable 



Case 5 is denned by conditions (|58j) and (JSHj). Using the phase j3 of the definition (J52j) . 
they read 

^— ^ + ^— J - + ±— >- = 0, (80) 



/i 2 



j"3 



"(^i) 2 + e 2 ^ (?7 Tl ) 2 " 










- + - 




- + - 



/il /i 2 /i3 

These two equations lead to 

fit _ ^2 (^2 + e^U r2 ) (if 2 , + e 2 ^) - U el (%t + e^U Tl ) (lf 2 2 + e 2 ^U, 



0. (81) 



Ms 

/'2 



f> e3 (^3 + eVUra) 




+ e 2i m 2 2 ) 


- Ue2 (U„2 + eWUrt) 




+ e 2i ?U 2 3 ) 


Ue2 (%2 + e^U r2 ) 


0% 




- Uei + e^U Tl ) 




+ e 2 ^U 2 2 ) 



^ U el + e#U Tl ) (U% + e 2 ^U 2 3 ) - U e3 (p^ + e^U r3 ) (if 2 , + e^U 2 Tl ) ' 



(82) 



33) 



In each of cases 6-9 there is a set of two conditions of the type {M. v l 



a/3 



and 



i pa 



0. By solving the system of these two equations one obtains 



to 
fH 



(U 12 U S 2 + e^U p2 il a2 ) U al U pi - (u^Usi + e^U pl U al ) U a2 U, 
(% 3 U S3 + e^U p3 U a3 ) U a2 Up 2 - (U l2 U S2 + e^U p2 U a ' 
(U J2 U S 



a2 J U a3 Up 3 

'62 + e^U p2 U a2 ) U al U pi - (tj^Usi + e*<>U pl U al ) U a2 Up 2 
(p^Un + e^U pl U al ) U a3 U 03 - (p l3 U S3 + e^U/aUrt) U al U t 



34) 
35) 



0i 



where ip = $ p + — $ 7 — $5. The unphysical phase ip arises because the relative phase 
in the equality (.M" 1 ) s = — (M.~ l ) is arbitrary due to the rephasing freedom. 



The right-hand sides of equations (J%2|) - (J%KJ) are functions of the observables 9\ 2 , 9 23 , 
9i3, and 5, plus an unphysical phase — f3 in equations (JS^j) and (J83|) . ip in equations (jSljl 
and 1(55)1 . 

We have used a numerical approach to cases 5-9, similar to the one in We have 
numerically generated random parameter sets of the four observables 6*12, 6*23, 9± 3 , and 5, of 
the unphysical phase, and of the mass m 3 . For each of these sets, the neutrino mass ratios 
mi/m 3 and m 2 /m 3 have been computed by taking the absolute values in the corresponding 
equations for fi\j fi-& and fi 2 / fi 3 ] by taking the absolute values of those equations we avoid 
considering the two observables on which we have hardly any experimental grip — the 
Majorana phases and Q. Whenever the obtained neutrino mass ratios, together with 
the inputted m 3 , were consistent with the experimental values for m\ — m\ and \m\ — m\\, 
we plotted the corresponding physical point in a scatter plot. The resulting scatter plots 
of the allowed regions for each model are presented in figures We have generated 
100,000 random parameter sets for each of the cases 5, 6, 7, and 9; for case 8 we have used 
2,000,000 random parameter sets because we have found that only a very small percentage 
of them turned out to be allowed. 
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We present in figures HHS1 for each of cases 5-9, respectively, two types of scatter plots 
of the allowed regions: in the sin 2 2# 23 -|[y e3 | plane, and in the m 2 -m3 plane. 

The neutrino mass spectrum in cases 5 and 7 may be of any possible type: normal, 
inverted, or quasi-degenerate. Case 6 typically displays a normal spectrum (m 3 > m 2 ). 
Cases 8 and 9 have inverted neutrino mass spectra (m 3 < m 2 ), as derived already in sec- 
tionElfrom the condition (A / l~ 1 ) ee = alone. Notice that case 8 has a distinct prediction: 
it requires the atmospheric neutrino mixing to be very close to maximal, and \U e ^\ is also 
usually close to its upper bound; this, too, we had already anticipated analytically in 
section 03 

For all cases 5-9, \U e ^\ cannot vanish, but it may be very small. We have performed 
a dedicated search of points with very small \U e ^\ by constructing random parameter 
sets with very small inputted #13. We thereby found the lower bounds on \U e ^\ given in 
tabled One sees that, in all cases but in case 8, \U e %\ can attain values so low as to be 





case 5 


case 6 


case 7 


case 8 


case 9 


\Uel\ > 


0.001 


10~ 6 


0.0005 


0.04 


0.0003 


(m) ee (eV) > 


0.0001 


10~ 5 


0.0001 


0.01 


0.01 



Table 1: Lower bounds of \U e s\ and on {m) ee , for each case. 



experimentally undistinguishable from zero. 

We have also investigated the decay rate for neutrinoless double-beta decay in each 
model. That rate is controlled by the effective Majorana mass 

(m) ee = mic 2 3 c 2 2 e ie + m 2 c 2 3 s 2 2 + m 3 s 2 3 e i(n ~ 2<5) . (86) 

The values of the Majorana phases are extracted from the equations for //1///3 and ^ 2 /^ 3 . 
We present scatter plots of (m) ee as a function of m 2 , for each case, in figures El and [7| 
One sees in those figures that (m) ee is at most 0.3 eV in cases 5-7; in cases 8 and 9 the 
upper bound is about an order of magnitude smaller. Lower bounds of (m) ee are given 
in tabled since cases 8 and 9 have an inverted neutrino mass spectrum, the lower bound 
on (w)ee is not so suppressed there as in cases 5-7. 



7 Summary 

The type-I seesaw mechanism has three sources of lepton mixing: the charged-lepton 
mass matrix Me, the neutrino Dirac mass matrix Md, and the mass matrix Mr of the 
right-handed neutrinos. It is interesting to choose Mr as the only source of mixing. A 
symmetry reason for this choice is provided by the assumption of conservation of the family 
lepton numbers L a (a = e, /1, r) in all terms of dimension 4 in the Lagrangian, because 
in that case M^ and Md are automatically diagonal; the Majorana mass terms of the 
right-handed neutrinos break the L a softly and we are therefore allowed to obtain lepton 
mixing from Mr [7j. This framework admits an arbitrary number of Higgs doublets, as 
their Yukawa-coupling matrices are all diagonal, hence flavour-changing Yukawa neutral 



12 



interactions among the charged leptons are strongly suppressed ^3]. in this framework, 
it looks sensible to consider relations among the matrix elements of the inverted neutrino 
mass matrix Ai'J 1 [10J. 

In this paper, we have used as a starting point a fi—r interchange antisymmetry in 
M.y 1 - We have shown that this can be obtained by imposing on the full Lagrangian a 
symmetry of the Z 4 type. However, such a /i-r antisymmetric matrix M.~ x is singular, 
therefore /i-r symmetric "perturbations" must be added to it. We have done so by 
introducing complex scalar SU(2) x £7(1) invariants with Yukawa couplings to the right- 
handed neutrinos. These scalar fields have non-trivial lepton numbers L a , and there 
are four basic possibilities to introduce such singlets in accord with the above-mentioned 
— see section 2. The VEVs of these scalar SU(2) x U(l) invariants generate the /i-r 
symmetric "perturbations" in M.~ x . 

Each of the four basic possibilities leads to a four-parameter neutrino mass matrix 
(cases 1, 2, 3, and 4 in section 2) which we have shown to be ruled out by the data. As 
a next step, we have considered the six combinations of the four basic possibilities. One 
of them leads to a five-parameter mass matrix (case 10) and is ruled out as well. On the 
other hand, the remaining cases 5, 6, 7, 8, and 9, which have neutrino mass matrices with 
six parameters, are physically viable, as we have shown by numerical analysis. We have 
refrained from combining more than two of the four basic possibilities, since the ensuing 
neutrino mass matrices presumably have negligible predictive power. 

The five viable cases have in common that they allow s\ 3 to be as large as the present 
upper bound. The numerical scan in the parameter space of cases 5, 6, and 7 shows 
a preference for a hierarchical spectrum, although quasi-degenerate or inverted spectra 
cannot be ruled out. The neutrino mass spectrum in cases 8 and 9 has inverted hierarchy; 
this follows solely from the condition (A^~ 1 ) ee = 0, which has been analytically analyzed in 
section 5, an analysis which is borne out by our numerical results — see figures 4 and 5. The 
most predictive neutrino mass matrix is the one of case 8, because it leads to practically 
maximal atmospheric mixing, correlated with a large sf 3 — see the analytical discussion of 
conditions (|57j) and (JBUj) in section 5, and the numerical result in figure 4. 

The research presented in this paper displays once again the usefulness of conditions on 
the matrix elements of the inverted neutrino mass matrix Ai' 1 . Indeed, in the framework 
employed here, it is most natural that symmetries imposed on the Lagrangian lead to 
relations in M.~ l . In this paper we have used /i-r interchange antisymmetry, instead of 
the well-known /i-r interchange symmetry, as the basic ingredient to obtain viable and 
predictive neutrino mass matrices. 
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Figure 1: Scatter plots of the allowed region in case 5, where 100,000 random parameter 
sets are generated. 
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Figure 2: Scatter plots of the allowed region in case 6, where 100,000 random parameter 
sets are generated. 
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Figure 3: Scatter plots of the allowed region in case 7, where 100,000 random parameter 
sets are generated. 
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Figure 4: Scatter plots of the allowed region in case 8, where 2,000,000 random parameter 
sets are generated. 
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Figure 5: Scatter plots of the allowed region in case 9, where 100,000 random parameter 
sets are generated. 




Figure 6: Scatter plots of the (m) ee of neutrinoless (5(5 decay as a function of m 2 for 
cases 5, 6, and 7. 
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Figure 7: Scatter plots of the (m) ee of neutrinoless (3(3 decay as a function of m 2 for 
cases 8 and 9. 
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